Spontaneously modulated spin textures in a dipolar spinor Bose-Einstein condensate 
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Helical spin textures in a 87 Rb F — 1 spinor Bose-Einstein condensate are found to decay spon- 
taneously toward a spatially modulated structure of spin domains. This evolution is ascribed to 
magnetic dipolar interactions that energetically favor the short-wavelength domains over the long- 
wavelength spin helix. This is confirmed by eliminating the dipolar interactions by a sequence of rf 
pulses and observing a suppression of the formation of the short-range domains. This study confirms 
the significance of magnetic dipole interactions in degenerate 87 Rb F = 1 spinor gases. 



In a wide range of materials, the competition between 
short- and long-range interactions leads to a rich land- 
scape of spatially modulated phases arising both in equi- 
librium and as instabilities in non-equilibrated systems 
[E 01- I n classically ordered systems such as magnetic 
thin films Q and ferrofluids 0, short-range ferromag- 
netic interactions are commonly frustrated by the long- 
range, anisotropic magnetic dipolar interaction, render- 
ing homogeneously magnetized systems intrinsically un- 
stable to various morphologies of magnetic domains [B[ . 
Long-range interactions are also key ingredients in many 
models of strongly correlated electronic systems @ and 
frustrated quantum magnets [7[. 

In light of their relevance in materials science, strong 
dipole interactions have been discussed as an important 
tool for studies of many-body physics using quantum 
gases of atoms and molecules, offering the means for 
quantum computation simulations of quantum mag- 
netism [§| and the realization of supersolid or crystalline 
quantum phases [Tol . However, in most ultracold 



atomic gases, the magnetic dipolar interaction is dwarfed 
by the contact (s-wave) interaction. Hence, experimental 
efforts to attain dipolar quantum gases have focused on 
non-alkali atoms, notably 52 Cr with its large magnetic 
moment 12 1, and on polar molecules [l3j ]. 

In this Letter, we demonstrate that magnetic dipole 
interactions play a critical role in the behaviour of a 
quantum degenerate F = 1 spinor Bose gas of 87 Rb. 
In this quantum fluid, s-wave collisions yield both a 
spin-inde pen dent and a spin-dependent contact interac- 
tion [HI ll5L Il6j |. with strengths proportional to a = 
(2a2 + ao)/3 and Aa — (02 — ao)/3, respectively, where 
the scattering length a f describes collisions between par- 
ticles of total spin F. In 87 Rb, with a Q (o 2 ) = 5.39 (5.31) 
nm, the spin-dependent contact interaction is far weaker 
than the spin-independent one; nevertheless, it is a criti- 
cal determinant of the ma gnet ic properties of degenerate 
F = 1 87 Rb gases [1_7|, UM, LL9j . The magnetic dipole in- 
teraction strength may be parameterized similarly by a 
length ad = fiog F [i B m/ (I2irh 2 ) , where /j,q is the perme- 
ability of vacuum, gp = 1/2 the gyromagnetic ratio, \ib 
the Bohr magneton and m the atomic mass [2(J . Given 
ad I Aa = 0.4, the F = 1 spinor Bose gas of 87 Rb is an 



essentially dipolar quantum fluid 21|. 

In our experiment, the influence of dipolar interactions 
on the spinor gas is evidenced by the spontaneous dis- 
solution of deliberately imposed long-wavelength helical 
spin textures, in favor of a finely modulated pattern of 
spin domains. We ascribe the emergence of this modu- 
lated phase to the magnetic dipole energy that disfavors 
the homogenously magnetized state and drives the fluid 
toward short- wavelength spin textures. To test this as- 
cription, we re-examine the behavior of spin helices in 
condensates in which the dipolar interaction is eliminated 
using a rapid sequence of rf pulses. The suppression of 
the modulated phase observed in this case confirms the 
crucial role of dipolar interactions. 

For this work, spin-polarized 87 Rb condensates of up 
to 2.3(1) x 10 6 atoms in the \F = 1, mp = —1) hyperfine 
state and at a kinetic temperature of T ~ 50 nK were 
confined in a single-beam optical dipole trap character- 
ized by trap frequencies (uj x ,oj y ,uj z ) — 2ir x (39,440,4.2) 
s _1 . The Thomas- Fermi condensate radius in the y (ver- 
tical) direction (r y = 1.8 /im) was less than the spin 
healing length £s — (87r Aa no)~ 1//2 = 2.4 fim where 
tiq = 2.3 x 10 14 cm~ 3 is the peak density of the con- 
densate. This results in a spinor gas that is effectively 
two-dimensional with regard to spin dynamics. 

The condensate was transversely magnetized by apply- 
ing a 7r/2 rf pulse in the presence of an ambient magnetic 
field of Bq — 165(5) mG aligned to the z axis. Stray mag- 
netic gradients (curvatures) were canceled to less than 
0.14 mG/cm (4.3 mG/cm 2 ). A helical spin texture was 
then prepared by applying a transient magnetic field gra- 
dient dB z /dz for a period r p = 5 - 8 ms. Larmor pre- 
cession of the atomic spins in this inhomogeneous field 
resulted in a spatial spin texture with a local dimension- 
less spin of F = cos(kz + uji£)x + sin(Kz + u;i,t)y, where 
k = (gp^B / h)(dB z / dz)r p z is the helix wavevector. The 
fast time variation at the 10^/2^ ~ 115 kHz Larmor pre- 
cession frequency will be henceforth ignored by consider- 
ing the spin at a particular instant in this rapid evolution. 
The helix pitch A = 27t/k ranged between 50 and 150 fim. 
Given A » ^s, the kinetic energy per atom in this spin 
texture, E K = H 2 K 2 /4m, was always negligible compared 
to the ferromagnetic contact interaction energy 
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FIG. 1: Spontaneous dissolution of helical textures in a quantum degenerate Rb spinor Bose gas. A transient magnetic 
field gradient is used to prepare transversely magnetized (b) uniform or (a, c) helical magnetization textures. The transverse 
magnetization column density after a variable time T of free evolution is shown in the imaged x — z plane, with orientation 
indicated by hue and amplitude by brightness (color wheel shown), (b) A uniform texture remains homogeneous for long 
evolution times, while (c) a helical texture with pitch A = 60 /tm dissolves over ~200 ms, evolving into a sharply spatially 
modulated texture. 



The helical spin texture was then allowed to evolve in a 
homogenous magnetic field for a variable time before the 
vector magnetization was measured using a sequence of 
non-destructive phase contrast images. Because of Lar- 
mor precession, a rapid sequence of images taken with 
circularly polarized light propagating along the y direc- 
tion can be analyzed to determine the column-integrated 
magnetization perpendicular to the ambient field 23l.l24j|. 
with vector components M XtV = (gF^B)nF XtV where h 
the column number density. Subsequent to this imaging 
sequence, a ir/2 pulse was applied to rotate the longitu- 
dinal spin F z into the transverse spin plane, and a second 
sequence of images was obtained. A least-squares algo- 
rithm comparing data from the two imaging sequences 
allowed the longitudinal magnetization M z to be deter- 
mined [25|. 

The evolution of helical spin textures is portrayed in 
Fig. [T] While uniform spin textures (A 3> 2r z ) remained 
homogenous for long times, helical textures (A < 2r z ) 
spontaneously develop short-wavelength modulations of 
the magnetization. This modulated phase is character- 
ized by spin domains with typical dimensions of ~ 10 /im, 
much smaller than the pitch of the imprinted helix, with 
the magnetization varying sharply between adjacent do- 
mains. The spatial modulations nucleated in regions that 
varied from shot to shot but gradually grew to encompass 
the entire condensate. 

To quantify this behaviour, we considered the power 
spectrum of the spatial Fourier transform of the vec- 
tor magnetization, \M.(k x ,k z )\ 2 , where (k x ,k z ) is the 
spatial wavevector in the image plane. This spectrum 
was found to consist of two distinct components: a cen- 
tral component that quantifies the long-range order of 
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FIG. 2: Power spectrum of the spatial Fourier transform and 
the two-point correlation function G(x, z) for the initial spin 
helix (a, c) and the spontaneously modulated phase (b, d). 
These data are derived from the same image sequence shown 
in Fig. 1 (d). The images (a, c) correspond to an evolution 
time T — ms while (b, d) correspond to an evolution time 
T — 250 ms. The short-range spatial order is defined as the 
integrated spectral power in the annular region shown in (b). 



the helical texture, and a second concentration of spec- 
tral power at a discrete set of wavevectors of magnitude 
femod — 27r/(10 /im) representing the short-range order of 
the final modulated texture. After subtracting out the 
background representing image noise, we divided spatial 
Fourier space into regions indicated in Fig. [5] and defined 
the integrated spectral power in the central region (annu- 
lar region) as the parametrization of long-range (short- 
range) spatial order in the quantum fluid. 
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FIG. 3: Growth of the spontaneously modulated phase (•) 
coincides with a reduction in the integrated energy in the low 
spatial frequency region (■). The data shown correspond to 
an initial helical pitch of 60 /im. Inset: The initial growth rate 
7 of the modulated phase as a function of the helix wavevec- 
tor. These were extracted from linear fits of the short-range 
order parameter at short evolution times. 



The formation of the spontaneously modulated tex- 
ture is reflected in the reduction of the long-range order 
parameter and the concomitant rise of the short-range 
order parameter (Fig. [3]). During this process, the total 
spectral power was found to be roughly constant indi- 
cating that the bulk of the quantum fluid remains fully 
magnetized even as the long-range order is reduced. The 
growth rate 7 of the short-range order parameter deter- 
mined from such data was found to rise monotonically 
with the wavevector k of the initial helical texture. While 
the long-range order was found to decrease after suffi- 
ciently long evolution times even in condensates prepared 
with nearly uniform magnetization, we note that stray 
magnetic field inhomogeneities of ~5 /iG across the axial 
length of the condensate would by themselves produce 
a helical winding across the condensate over a period of 
300 ms, constraining our ability to test the stability of 
homogenous spin textures. 

Another measure of the spontaneous short-range mod- 
ulation in the condensate is the appearance of polar-core 
spin vortices throughout the gas. Such vortices were 
identified as in Ref. [l^ | by a net winding of the trans- 
verse magnetization along a closed two-dimensional path 
of non-zero magnetization in the imaged gas. The num- 
ber of identified spin vortices was roughly proportional to 
the short-range order parameter, with no vortices iden- 
tified in the initially prepared spin helix and up to 6 
vortices/image identified in the strongly modulated tex- 
ture produced after free evolution. In each instance, the 
number of vortices with positive and negative charge was 
found to be approximately equal. 

A striking feature in the evolution of spin textures is 



the significant rise in the kinetic energy of the condensed 
atoms, reaching a value of fi 2 ^ od /8m = h x 6 Hz per 
atom given that roughly half the spectral weight of the 
texture's magnetization is at the wavevector fc mo d- One 
expects the total energy per atom in the condensate to 
be constant during this evolution, or even to diminish 
through the transfer of energy to the non-condensed por- 
tion of the gas. Yet, in examining the energy of the 
initially prepared spin helix, we find the local contact- 
interaction energy is minimized, the quadratic Zeeman 
energy is just q/2 — h x 1 Hz at the ambient mag- 
netic field, and the kinetic energy of the spin helix is 
just E K /2 < h x 0.5 Hz for a helix pitch of A > 50 /im. 

This apparent energetic deficit of the spin helix state 
can be accounted for by the magnetic dipole interac- 
tion. The on-axis magnetic field produced by a spin 
helix in an infinite axial column of gas with a gaus- 
sian transverse density profile can be simply calculated. 
From this calculation, we estimate that a gas with uni- 
form transverse magnetization possesses an excess of 
Ed = Mo fl^ Ms n o/ 2 ~/ix5Hz compared to the energy of 
a tightly wound helix, a figure that closely matches the 
excess kinetic energy of the finely modulated texture. 

To confirm the role of magnetic dipolar interactions 
in the evolution of these spin textures, we employed a 
modification of the NMR technique of spin-flip narrow- 
ing (2(| to eliminate effectively the dipolar interactions. 
The interaction energy of two magnetic dipoles sepa- 
rated by the displacement vector r is proportional to 
Fi -F 2 — 3(f-Fi)(f -F 2 ). If both dipoles experience rapid, 
common rotations that evenly sample the entire SO(3) 
group of rotations, the interaction energy will average 
to zero regardless of the relative orientations of the spin 
vectors Fi^ and of the displacement vector r. We note 
that such spin rotations also annul the quadratic Zeeman 
shift. However, since the character of the spontaneously 
modulated phase was observed to be unchanged as q was 
varied over a factor of 5 (0.4 < ^ q/2 < 2 Hz), we ignore 
this small difference. 

Experimentally, after the initial spin texture had been 
prepared as before, we effected such spin rotations by ap- 
plying a rapid sequence of 7r/2 rf pulses to the Larmor 
precessing atoms at random intervals, and, thus, along 
random rotation axes, at a mean rate of I - 2 kHz. The 
spin helix was allowed to evolve under the constant ac- 
tion of these dipole-cancelation pulses until the pulses 
ceased and the sample was imaged as described earlier. 
As shown in Fig. [4j in the absence of dipole interactions 
a spin helix was still prone to spontaneous spatial modu- 
lation, consistent with the predictions of dynamic insta- 
bilities made in Refs. 27 1 and [28j]. However, the depth 
of this modulation was significantly suppressed, as seen 
from the smaller spectral weight of short-range magnetic 
order. This suppression was less evident for large he- 
lix wavevectors (n > 110 mrad//im), presumably due to 
the faster growth of the spontaneously modulated phase. 
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particular interest is the influence of such anisotropic in- 
teractions on the true ground state of these spinor gases 
and the question of whether the self-organized modulated 
spin texture represents an unforeseen equilibrium phase 
of this dipolar quantum fluid. 
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C. H. Greene and A. Lamacraft. This work was sup- 
ported by the NSF, the David and Lucile Packard Foun- 
dation, and DARPA's OLE program. Partial personnel 
and equipment support was provided by the Division of 
Materials Sciences and Engineering, Office of Basic En- 
ergy Sciences. S. R. L. acknowledges support from the 
NSERC. 



FIG. 4: Spatial power spectra (a) without or (b) with the 
application of rapid rf pulses during the 200 ms evolution 
following the preparation of the helical spin texture, are av- 
eraged over four experimental repetitions. Eliminating the 
dipolar interactions suppresses the short-range spatial mod- 
ulation of the final spin texture, (c) The average number of 
detected spin vortices and (d) the short-range order param- 
eter are shown vs. evolution time. Square symbols indicate 
data obtained when the dipolar interaction is spatially aver- 
In all cases, the initial helix pitch was A = 80 /an. 



The reduction of the excess kinetic energy of the final 
spin texture in the spin-rotation-averaged sample sup- 
ports our identification of the dipole energy as its source. 

Finally, we note the distinct six-fold structure in the 
spatial Fourier spectrum of the spontaneously modulated 
phase (Fig.^b)). The interpretation of this structure is 
aided by considering the spatial correlation function of 
the magnetization, which we define as 



G(Sr) 



J2 r M(r + Sr) ■ M(r) 
(9fHb) 2 E r fl ( r + <5r)n(r) 



(1) 



For the initial spin helix, this correlation function shows 
the long-range sinusoidal variation of the transverse mag- 
netization at the corresponding pitch. In contrast, the 
final texture shows significant short-range correlations in 
both the x and z directions, with regions of opposite 
alignment arranged in the form of a checkerboard with a 
lattice spacing l m ~ 10 ^m (Fig. 2(d)). While these cor- 
relations are strongest at short range, they persist, albeit 
with diminished strength, even for separations Sr l m . 
This lattice structure is suppressed in helices evolving 
under the active cancelation of dipolar interactions. 

In conclusion, we observe a magnetic-dipole-mediated 
instability resulting in the emergence of spontaneously 
ordered spin domains in a quantum degenerate spinor 
Bose gas. The demonstration of the significance of mag- 
netic dipolar interactions in a 87 Rb spinor gas presents 
a new arena for the study of dipolar quantum gases. Of 
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